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Abstract. We give a conjecture for the expected value of the optimal k- 
assignment in an m X n-matrix, where the entries are all exp(l)-distributed 
random variables or zeros. We prove this conjecture in the case there is a zero- 
cost k — 1-assignment. Assuming our conjecture, we determine some limits, 
as k = m = n — » oo, of the expected cost of an optimal ra-assignment in 
an n by n matrix with zeros in some region. If we take the region outside a 
quarter-circle inscribed in the square matrix, this limit is thus conjectured to 
be tt 2 /24. We give a computer-generated verification of a conjecture of Parisi 
for k = m = n = 7 and of a conjecture of Coppersmith and Sorkin for k < 5. 
We have used the same computer program to verify this conjecture also for 
k = 6. 



1. Introduction 

Suppose we are given an m by n array of nonnegative real numbers. A k- 
assignment is a set of k entries, no two of which are in the same row or the same 
column (we "assign" rows to columns, or vice versa). The cost of the assignment 
is the sum of the entries. A /c-assignment is optimal if its cost is minimal among 
all fc-assignments. 

We let Fk(m,n) denote the expected cost of the optimal fc-assignment in an m 
by n array of independent exponentially distributed random variables with mean 
1. 

Aldous [A95] has shown the existence of lim„^oo F n (n, n) — c. According to 
a conjecture of Mezard and Parisi [MP85|, c = ir 2 /6. Moreover, Parisi |P9§[ ] has 
conjectured that F n (n,n) = l + l/4+l/9 + -- - + l/n 2 . This conjecture has been 



verified for n < 7. For n < 6 this can be found in [ AS99|, and the n — 7 case follows 
from the calculations in the appendix to this paper. 
In [ CS98| , a more precise conjecture was formulated. 



Conjecture 1.1. 



F k (m,n) 



E 

i,j>0 
i+j<k 



1 



(m-i)(n- j) 



It is also shown in CS98 ] that this conjecture is in accordance with the conjecture 
of Parisi if k = m = n. We have verified Conjecture 1.1 for k < 6, (and arbitrary m 
and n). The cases k = 3, 4 are treated in | AS99] and k — 5 is given in the appendix. 
For k = 6, about 10,000 cases were needed. 
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1.1. The Main Conjecture. As is shown in ]CS98| , [ |AS99| , and in Section [| 
the calculation of Fk(m, n) can in some (but not all!) cases be done recursively by 
reducing the problem to several assignment problems for matrices which have zeros 
in certain positions. 

Let Z be a finite set of matrix-positions, that is, 

Z = • • ■ , (i r ,3r)}, 

for some positive integers ,i r ,jr- Suppose m > max(ii, ...,i r ) and n > 

max(ji, ...,j r ). We define Fk,z{Tn,n) to be the expected cost of the optimal k- 
assignment in an to by n matrix with zeros in the positions belonging to Z, and 
independent exponential random variables with mean 1 in all other positions. 

We will consider sets of rows and columns in the m x n-matrix. A set a of 
rows and columns is said to cover Z if every entry in Z is on either a row or on 
a column in a. A covering with s rows and columns will be called an s-covering. 
k — 1-coverings will be of particular importance. By a partial k — 1-covering of Z, 
we mean a set of rows and columns which is a subset of a k — 1-covering of Z . Let 
Pi t j be the probability that a set of i rows and j columns, chosen from uniform 
distribution on all such sets, is a partial k — 1-covering. 

We now present a conjecture for the values of Fk i z(m, n). 



Conjecture 1.2 (Main Conjecture). 

F kiZ (m,n) = 



E 

i,3>0 
i+j<k 



Pi,j 



(to — i)(n — j) 



(1) 



Example 1.3. When Z is empty, Conjecture 1.2 clearly specializes to Conjecture 







Example 1.4. When k = 2 and Z has just one element, then pi^o = po.i — 
1/n and pq.q = 1, so according to Conjecture 1.2 we get 

1 1/to 1/n -1 1 1 

\ ! 1 '- = 1 1 

ran (to — \)n m(n — 1) to 



F2,z(m,n) 



(m — V)n m(n — 1) 



which is also the correct value, see Example 2.4 below. 

In fact, we prove in Section [3] that Conjecture 1.2 holds under the assumption 
that Z contains a k — 1-assigmnent, which is our main theorem. 

Let P be the poset of all intersections of k — 1-coverings of Z, ordered by reverse 
inclusion, and with an artificial minimal element corresponding to the empty 
intersection. We prove in Section |^ that the following is an equivalent form of the 
Conjecture 1.2: 



Fk,z{m,n)= 



*sp\o 



(m - i a )(n - j a ) 



(2) 



where fj-p(a) is the Mobius function of the interval (0,a) in P and i a ,j a are the 
number of rows and columns in a respectively. 

In Section || we describe our tools to recursively compute Fk,z(m, n), which we 
use both for the theoretic proofs and for the algorithmic computational results in 
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the Appendix. In Section [| we discuss some consequences of Conjecture 1.2. One 
such consequence is that we always can write 

F k ,z(m,n) = ^ — 



(m — i)(n — j) 
i+ 3 <k y n ■>> 



for some integers 6jj independent of m and n. We also prove that Conjecture 1.2 



implies a conjecture by Olin [092 1 on the probability that the smallest element in a 
row is used in the optimal assignment. In Section ^ we prove that some asymptotic 



results follows from Conjecture 1.2 



Acknowledgement: We thank Mireille Bousquet-Melou and Gilles Schaeffer for 
introducing us to the problem. 

2. Computation of F k> z (m, n) 

If A is a (random) m by n matrix, and u is an assignment, then we let ^ U (A) 
be the cost of u, that is, 

If k is an integer not greater than m or n, we let ^ k (A) be the cost of the optimal 
fc-assignment in A. 

We say that a set a of rows and columns is an optimal covering of a set Z, if a 
covers Z, and a has minimal cardinality among all coverings of Z. Our calculations 
of Fk t z{m, n) are based on the following theorem. 

Theorem 2.1. Let k < m,n be positive integers. Let A be an m by n matrix, with 
zeros in a certain set Z of positions, and (possibly random) positive values outside 
Z . Let a be an optimal covering of Z . Let e be a positive real (possibly random) 
number, not greater than any of the entries not covered by a. Let A' be the matrix 
obtained from A by subtracting e from every entry not covered by a, and adding e 
to every doubly covered entry of A. Then 

* k (A) = e-(k-\a\) + * k (A'). 

To prove this theorem, we need the following lemma, which we prove in Section 



2.1 



Lemma 2.2. If a is an optimal covering of a subset of Z , then there is an optimal 
k-assignment which intersects every row and every column of a. 



Proof of Theorem 2A_. Let u be an optimal fc-assignment of A' intersec ting every 



row and column of a. Such an optimal assignment exists by Lemma Ifl. Note 
that subtracting e from every entry not covered by a, and adding e to every doubly 
covered entry is the same thing as first subtracting e from all entries not covered 
by the rows of a, and then adding e to all the entries covered by the columns of a. 
Let i and j be the number of rows and columns in a, respectively. Then 

*u(4) =V u (A')+e-(k-i)-e-j = 

= e ■ (k - H) + *„(A') = e ■ (k - H) + 9 k (A'). (3) 

For every fc-assignment v we have 

*v(A) >e-(k- H) + 9 V (A') >e-(k- \a\) + $ k (A'). 

Hence u is an optimal fc-assignment also in A, and ^ k (A) is given by (3). □ 
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In the typical use of the theorem, e is the minimum of all the non-covered entries, 
which gives a new zero when we subtract e. The special case of this theorem where 
there are no doubly covered entries or all such entries are known not to be in 
the optimal assignment, was proved and used in [CS9S] and [ AS99| to prove 



Conjecture 1.1 for k < 4 



Before proving Lemma 2.2, we give an example of how Theorem |2.l| is used. 
First we need to state some well-known properties of the exponential distribution. 
We say that a random variable is ea;p(7)-distributed if it is exponentially distributed 
with mean I. 

Lemma 2.3. Let X\, . . . ,X t be independent exponential random variables, with 
mean cti, ci2, • • • , a>t respectively. Let Y = min{X l7 . . . , X t }. Then 
1. Y is exponentially distributed with mean 

1 



l/oi + • • • + 1/at 

2. ifY = Xi, then Xj — Y is again exp(aj)- distributed, for i =^ j , and {Xj — Y : 
j =/= i} are independent, and independent ofY. 

3. if I is the random variable defined byY — Xi, then L is uniquely defined with 
probability 1 and 

l/a l 



Pr{I 



l/oi + --- + l/o t ' 



Example 2.4. Let us show how this is used to compute F2(m,n). First we take e 
to be the minimum of all elements of A, which by symmetry we can assume to be 
ai t \. We now subtract a^i from all entries in A to obtain the matrix A' with a zero 
in position (1, 1) and exp(l)-distributed random variables elsewhere. The covering 



a is here empty and e is exp(l/mn)-distributed. From Theorem 2.1 we get 



2 

F 2 (m,n) = VF 2Z ,{m,n), 

mn 

where Z' = {(1,1)}. 

To compute F 2 ^z'{jn, n), we have to find an optimal covering of Z' and we can 
choose a' — {n.}, the first row, see Figure [j]. We again take the minimum of the 
entries th at a re not covered by a', which means that e' is exp(l/(m — l)n). From 
Theorem 2.1 we get 



1 1 , . n — 1 

F2,z'(m,n) = — + -F 2Z »[m,n) H F 2 ,z»'Vn,n), 

(m — l)n n n 

where Z" = {(1, 1), (2, 1)} and Z'" = {(1, 1), (2, 2)}. The minimum is equally likely 
to occur in any of the entries that is not in the first row. Here Z" represents the case 
that the minimum occurs somewhere in the first column, while Z'" represents the 
case that the minimum occurs in another column. Hence we get only two different 
non-isomorphic cases. 

Fi,z in (™, n) = because Z'" contains a 2-assignment of zeros. Finally we use 



Theorem 2.1 on Z" . The only optimal covering is the first column, a" = {ci}. The 
minimum of the other entries is exp(l/m(rj — 1)), which we use as e". Note that 
a new zero in another column than the first will produce a 2-assignment together 
with one of the zeros in Z" . We get 

F 2 . Z "(m,n) = 1 

m(n — 1) 
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Substituting back we get 

F2,z'(m,n) 



as in Example 1.4, and 



mn (to — l)n m(n — 1) ' 



F 2 ,z{m 1 n) 



in accordance with Conjecture 1.2 




to 7i (m — l)n TO,(n — 1) ' 



Figure 1. The cases needed when computing F2(m,n). The 
shaded row and column are the coverings. 



Example 2.5. Let us also give an example where there is a doubly covered entry. 
Let k = 5 and Z — {(1, 2), (1, 3), (2, 1), (3, 1)}. This is the simplest case where the 
recursion will result in a matrix with an entry which is neither exp(l) nor zero. 

The only optimal covering a = {ri, ci} consists of row one and column one. As 
usual we let e be the minimum of the non-covered entries, which here is exp(l/(m — 
l)(n- 1)). We get 

3 1 
F 5 ,z (m, n) = — + — ^ J! *5 ( A *,3 ) ' 

(ij')e[2,ro]x[2,n] 



(to — l)(n — 1) (to — l)(n — 1) 



where A i ■ has zeros in Z U {(i, j)} and an entry which is the sum of two r.v.'s, 
a\ t x = cxp(l) + exp(l/(TO — l)(n — 1)). There are actually four different cases here, 
and as is shown in the appendix it is possible to continue the recursion in each case 
to compute z(m, n ) and finally F^(m, n). 







Figure 2. The smallest "problem case". The shaded row and 
column is the covering. 
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Note that in Section 2.3 we show a way to avoid getting sums of random variables 



as entries in this particular case. We can thus reduce the number of "problem 



cases", that is cases when the recursion in Theorem 2.1 gives us entries different 



from exp(l) and zero. But we have not found any way to avoid them completely. 



2.1. Proof of Lemma 2.2. We begin by citing a famous theorem of D. Konig: 

Theorem 2.6. If a set of matrix positions cannot be covered by k — 1 rows and 
columns, then it contains a k-assignment. □ 

Lemma 2.7. Let k < m,n be positive integers. Let A be an m by n matrix, with 
zeros in a certain set Z of positions, and (possibly random) positive values outside 
Z. Suppose that Z does not contain a k + l-assignment. Let a be an optimal 
covering of Z. Then every row and every column of a contains an element of every 
optimal k-assignment of A. 

Proof. Suppose that a contains rows 1, . . . , r and no other rows. Since a is optimal 
it follows by Konig's theorem that there is an r-assignment v in Z containing no 
element from the columns in a. Suppose (for a contradiction) that there is an 
optimal k- assignment u which does not use row 1. 

We now construct a sequence of matrix positions (all with zeros) as follows: Let 
v be the element of v which is in the first row. Suppose that we have defined 
Vq, . . . ,Vh and u\, . . . , Uh- Then if there is an element of u D Z in the same column 
as Vh, let this element be Uh+i, and let vt+i be the element of v which is in the 
same row as Uh+i- Since u does not contain any element from the first row, the 
sequences vo,vi,V2, ■ ■ ■ and ui, U2, ■ ■ • cannot contain any repetitions of the same 
element. Hence the sequence must end with an element Vh such that there is no 
element of u PI Z in the same column. 

We now consider two cases. Suppose first that no element in the column of 
Vh belongs to u. Since Z does not contain a k + l-assignment, the cardinality 
of a is at most k. Each row and column of a covers at most one element of 
u, and row 1 does not cover any element of u. Consequently there is an ele- 
ment (i,j) of u which is not covered by a, hence does not belong to Z. Then 
u \ {u%, . . . , Uh-i, (i, j)} U {vq, . . . , Vh) is a fc-assignment of smaller cost than u, a 
contradiction. If on the other hand there is an clement of u in the first 

column, then u \ {ui, . . . , it/i-i, («', j')} U {vq, . . . , Vh} is a fc-assignment of smaller 
cost than u. This contradiction proves the lemma. □ 



Proof of Lemma \2.\ This follows by a continuity argument. Let the values in the 
positions of Z not covered by a be 6, and let 5 tend to zero. Since there are only 
finitely many fc-assignments, there has to be an assignment which is optimal for all 



sufficiently small 5, which by Lemma 2.7 has to intersect every row and column of 



a. This assignment is also optimal for 5 = 0. □ 

Note that it is not true that every optimal fc-assignment has to intersect the rows 
and columns of an optimal covering of a subset of Z. 

2.2. Superfluous matrix elements. It can be of great help to know that a certain 
entry of the matrix cannot occur in any optimal assignment. Sufficient conditions 



to draw this conclusion were given in Lemma 26 of [ CS98 and Lemma 9 of [ AS99| 



The following lemma gives a necessary and sufficient condition. The condition here 



is in fact equivalent to the condition given in JAS99 
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Theorem 2.8. Let k < m,n be positive integers. Let A — (a^j) be an m by n 
matrix. Suppose that a set Z of entries are zero, and the remaining entries are 
random variables that can take arbitrary positive values. Suppose that £ Z. If 
every set ofk—1 rows and columns covering Z also covers then cannot 

belong to an optimal k- assignment. Conversely, if there is a k — 1- covering of Z 
which does not cover (i,j), then it is possible to assign positive values to the matrix 
entries not in Z in such a way that belongs to a unique optimal k-assignment. 

Proof. By the theorem of Konig, a set contains a fc-assignment if and only if it 
cannot be covered by fewer than k rows and columns. The minimal cost of a 
fc-assignment is of course the same as the minimal cost of a set containing a fc- 
assignment. In trying to minimize the cost of such a set, it suffices to consider sets 
containing Z, since the entries in Z do not increase the cost. 

Suppose that (i,j) is covered whenever Z is covered with fewer than fc rows and 
columns. Then from any set containing Z, which cannot be covered with fewer than 
fc rows and columns, we can delete and obtain another such set of smaller cost. 
This shows that can never belong to an optimal fc-assignment. 

Suppose on the other hand that we can cover Z with a set a of fc — 1 rows 
and columns, without covering Then we assign the value 1 to ajj, and to 

the entries not in Z which are covered by a, and values larger than fc to all other 
entries not in Z. Since m and n are not smaller than fc, the set of entries which are 
covered by a cannot be covered in any other way by fc — 1 rows and columns. Hence 
the set of entries which are < 1 cannot be covered with fc — 1 rows and columns. By 
Konig's theorem there is a fc-assignment of cost at most fc. However, since every fc- 
assignment must use at least one entry not covered by a, every fc-assignment of cost 
at most fc must use djj. Hence aij belongs to the unique optimal fc-assignment. □ 



In this situation, the fact that Fk.z{rn,n) = -Ffc.zu{(i,j)} n ) can sometimes 
simplify computation. It is possible that several elements have the property of 
in the above theorem. Indeed, as the following theorem shows, it is sometimes 
possible to dispose of an entire row or column. If r is a row, we let Z\r denote the 
set where all positions of Z in row r has been removed and similarly for a column. 

Theorem 2.9. Suppose that a row r is used in every covering of Z by k — 1 rows 
and columns (in particular, this is the case if at least fc elements in r belong to Z). 
Let Z' = Z\r. Then 

F kiZ {m,n)=F k _x,z>{m-l,n). (4) 

Similarly, if a column c belongs to every covering of Z with fewer than k rows and 
columns, and Z" = Z\c, then Fk,z{"rn,n) = Fk t z"(m,ri — 1). 

Proof. If r is used in every covering of Z by fc — 1 rows and columns, then by 



Theorem 2.8, no non-zero element in row r, no matter how small, is used in any 
optimal fc-assignment. Hence we do not change the cost of the optimal fc-assignment 
if we replace all non-zero entries in row r by zeros. Let A' be the matrix obtained 
from A by deleting the row r. Every fc-assignment of A must of course contain (at 
least) a fc — 1-assignment of A' . It is not a problem if Z contains an entry in row 
m, since we can always permute the rows in Z so row m becomes empty instead of 
row r. 



s 
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Conversely, if row r consists of only zeros, then every k — 1-assignment in A! can 
be extended to a fc-assignment in A of the same cost. Hence ^k(A) — ^k-i(A'). 
Passing to expected values, we obtain (^). 

The corresponding statement for columns is of course equivalent. □ 

2.3. The probability that a certain element belongs to the optimal assign- 
ment. In this section we prove a theorem relating the expected values of certain 
assignments to the probability that a certain element belongs to such an assign- 
ment. We show how this implies a second recursion which can be used to avoid 
some "problem cases" in the computation of Ffc^(m,n). 

Theorem 2.10. Let k < m,n be positive integers. Let A = (ajj) be an m x 
n-matrix, where a set Z of the entries are zero, and the others are independent 
exponentially distributed random variables with mean 1. Suppose that ^ Z, 

and let Z' = Z U {{i,j)}- Then the probability that (i,j) belongs to an optimal 
k-assignment in A is F^^zijn^n) — Fk^z'ifn^n). 

Proof. We condition on the values of all matrix elements except aij. Let X be 
the value of the minimal fc-assignment in A which does not use ajj. Let Y be the 
value of the minimal k — 1-assignment in A which does not use row i or column j. 
The optimal fc-assignment in A either contains and has value y + a^j, or does 
not contain and has value X. Hence the probability that a, ( j belongs to the 
optimal fc-assignment in A is equal to the probability that a^.j < X — Y . 
We wish to show that this is equal to 

Fk,z(m,n) - Fk,z'(m,n) = E(max(0, min(X - Y, cn t j)). 

If X < Y, then both are zero. If X > Y, then we let 5 = X - Y. Then 
the probability that a iy j is used in the optimal fc-assignment in A is 1 — e~ 5 . We 
compare this to 

r s 

E(min(<5, a tJ )) = S e - A + / te' 1 dt = Ser 5 + 1 - (S + l)e~ A = 1 - e~ s . 
Jo 

This proves the theorem. □ 



From this we can get another generalization of the basic recursion used in JCS98| 



and [AS99|, which is of some help when all the doubly covered entries are nonzero. 



However, none of the results in this paper depends on this recursion. 

Corollary 2.11. Let fc < m, n be integers. Let Z c [1, m] x [1, n], be a set of zeros. 
Let a be an optimal covering of Z and let R be the non-covered entries. Also let 
the random variable S be the number of positions not covered by a that are in the 
optimal assignment. Then 

E[S] J_ x ^ 

W + W\ 



(i,j)ER 

Proof. By linearity of the expected value we have 

E[S}= ^2 (Fk,z(m,n) - F kyZu ^^(m,n). 

□ 
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Example 2.12. Let us give an example of how Corollary [2. 11 can be used to 



[2.5| we 



simplify the calculations of F k-Z (m, ri). As in Example [2.5| we take fc = 5 and 
Z = {(1,2), (1,3), (2, 1), (3, 1)}. The optimal covering a is the first row and the 
first column. The expected number of non-covered elements used in the optimal 
k-assignment is 3+ the probability that a\.\ is in the optimal assignment which is 
F k ,z(m,n) - F kiZ u{(i,i)}(™>,ri). We get 

„ , N 3 + F kiZ (m,n)-F kzU {(ii\}(m,ri) 

F k ,z{m,n) = — + 

(m — l)(n — 1) 

+ 7 TT7 77 Fk,zu{(Lj)}(m,n). (5) 

(m— l)(n— 1 ), , , , 

(ij)£[2,m]x[2,n] 

This can then be rewritten as 

n , , 3 - F fc,ZU{(l,l)}(m,^) + E(i,j)e[2,mlx[2,n] F fc,^U{(i,^}( m ' n ) 

i< k ,z{m,n) = / ttt tt , 

(m — l)(n — 1) — 1 

which enables us to stay within the set of matrices with all entries exp(l) or zero. 
However, this recursion does not help in the case 

Z = {(1,1), (1,2), (1,3), (2,1), (3,1)}. 

For this we would need to know the probability that the zero in position (1,1) is 
in the optimal k-assignment. We know of no simple way to calculate this, and 
in general it might not even be well-defined, since in some cases with positive 
probability there are several optimal /c-assignments using the same non-zero entries, 
but different zeros. 

3. Proof of the Main Theorem 



In this section we give a proof of Conjecture 1.2 under the assumption that Z 
contains a k — 1-assignment. 

Let k, m, n and Z be as usual. Let B(m, n) be the Boolean algebra of all sets of 
rows and columns in A, ordered by reverse inclusion. We define a function / = f k ^z 
on B(m, n) by 



/(«) = 

Wc define g = g k ,z by 



1, if a is a partial k — 1-covering 
0, otherwise 



E 9ifi) = /(«)• 

0<a 



Lemma 3.1. Let P = P kt z{m,n) be the sub-poset of B{m,n) consisting of all 
intersections of k — \-coverings, together with an artificial minimal element. Then 



9 (a) 



-/if. (a), ifaeP 
0, ifai P 



Proof. Suppose that a is minimal contradicting the statement. Then a ^ P and 
g(a) ^ 0. Let a' be the intersection of all atoms of P which are smaller than a, 
that is, which are supersets of a. Then a' 6 P, and for every (3 € P, if ft < a, 
then ft < a'. Hence ^Zp <a g{ft) = 0, and consequently, g(a) — 0. This is a 
contradiction. □ 
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For a S B(m,n), we let i a and j a be the number of rows and columns in a, 
respectively. 

As we will prove, a reformulation of our main conjecture is the following: 



Conjecture 3.2. 



F k ,z(m,n) = ^2 



a£.B(m,n) 



(m - i a )(n - j a ) ' 



(6) 



Example 3.3. Let k = 3 and Z = {(1,1)}. We will compute ^3,2(771,71) using 
Conjecture 3.2 and Lemma 3.1. The k — 1-coverings of Z are all sets of two rows 
and columns including row one or column one. The poset P = P3,z(rn,n) and the 
Mobius function is shown in Figure ||. The sets with an index i, e.g. rirj, represents 
all m — 1 sets with 2 < i < m and similarly the sets with an index j represents all 
n — 1 sets with 2 < j < n. We get 



Fz,z{m, n) = (m - 1)- 



1 



+ (n-l) 



1 



1 



(777 — 2)77 

+ (™- 1 )t TT7 7T + ( n ~ 1 ) ' 



(777 — 1)(77 — 1) 

+ (m-l) 7 ^— + (77-1) ' 



(m-l)(n-l) (to- l)(n- 1) 
—m — 77 + 2 

+ 



777(77 — 2) (m — 1)ti 

— 777 — 77 + 2 777 + 77 — 2 



(777 — l)n 
1 



777(77 — 1) 777(77 — 1) 

1 1 



1 



(777 — 2)77 777(?7 — 2) (777 — 1)(77 — 1) 777(77 — 1) (777 — 1)77 77777 

which is the correct answer. 



— , (7) 



-m-n+2 




FIGURE 3. The poset P3 ) z(777,n) in Example |3 
function is written in boldface. 



The Mobius 



We now state a lemma, which will also be useful later. It should be known, but 
for completeness we give a proof. 
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Lemma 3.4. If m and r are positive integers, then 

(-l) r 

i=0 

Proof. We have 



m - i m( m r 1 ) 



m — I 1 — ' /n \ I 

i=0 i=0 u v 



■ly f x m -'- r J2(-i) 

Jo 



dx 



r 

r — i 



i=0 

l 

m — 1 — r 



(-l) r / x m -^- r {l-x) r dx. (8) 



This integral can be interpreted as a probability. Let Y\ , . . . , Y" m be independent 
random variables uniformly distributed in [0, 1]. Then the integral in (8) represents 
the probability that Yi < Y m -r for i = 1, ...,m — r — 1, and Yj > Y m - r for 
j = m — r + 1, . . . , to. If we sort the to numbers Yi, . . . , Y m according to size, every 
permutation occurs with equal probability 1/to!. The number of permutations such 
that Yi, i = 1, . . . , to — r — 1 occur before y m _ r , and Yj, j = m — r + 1, . . . , to all 
occur after Y m - r , is r\{m — r— 1)!. Therefore, the probability that Yi < Y m _ r < Yj 
for alli = l,...,m — r — 1 and j = to — r + 1, . . . , to, is 

H(to — r — 1)! 1 

ml = mC"- 1 ) ' 

This proves the lemma. □ 



Proposition 3.5. Conjecture l.i Conjecture 3.L 

Proof. Let C be the set of partial k — 1-coverings. For (3 E C, we define 



/Ma) = < 

Obviously, / = J2pec hi- Moreover 



1, if a = [3 
0, otherwise 



/3ec 

where c/^ is defined by 

^5/3(7) = //3(a)- 

We have 17/3(0/) = /i(/3, o/), the Mobius function of the interval (/3, a) in the 
Boolean algebra over the set of rows and columns (ordered by reverse inclusion). 
This interval is itself a Boolean algebra. Hence 

f(-l)l/»l-M, if a C{3 
I 0, otherwise 



9 pia) 
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For the right hand side of (0), we have 



g(a) ^ x ^ (_l)l/»l-M 



^— ' (m — i n )(n — (to — irv)(n— 7V,) 

V( 

3=0 




[by Lemma 3.4 



1 



f-,TO( m - 1 W»r 1 ) 



The factor 



( m )( n ) 



yip' yon* 

can be interpreted as the probability of (3 under the uniform probability measure 
on all sets of ip rows and jp columns. Hence (9) equals ([!]). □ 

We need the following lemma. 

Lemma 3.6. Suppose Z contains a k — l-assignment u. Let Z 1 = Z U {(i,j)}> 
where (i,j) lies in a column not used by u. Then 



9k,z>{a) 



gk,z(ct), if a contains row i 
0, otherwise 



Proof. It follows from the assumptions that in order to cover Z' with k — 1 rows 
and columns, one must use row i. This implies that Pk^z'im, n) consists of precisely 
those elements of Pk,z{rn, n) which contain row i. Since this means that Pk,z' (jn, n) 
is an order- ideal in Pk,zijn, n), the Mobius function on Pk,z' ( m , n) is the restriction 
of the Mobius function on Pk.z(ni,n) to Pk,z' {m> n). Now the statement follows 



from Lemma 3.1. □ 



We now prove the formula of Conjecture 3.2 under the condition that Z contains 
a k — l-assignment. 

Theorem 3.7. If Z contains a k — l-assignment, then 



F k , z (m,n)= £ ^ -. (10) 

Proof. If Z contains a fc-assignment, the statement is true. We will therefore assume 
that Z does not contain a k- assignment. We will use induction on k, and for 
a specific value of k, we will also use induction on the number of elements not 
belonging to Z. Hence we will assume that the statement has been proved for 
Fk,z', for every superset Z' of Z. 

Notice that if there is a row r which belongs to every k — 1-covering of Z, then 
by Theorem [2~<i| , 

Fk,z(m,n) = F k _ hZ \ r (m - l,n). 
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By induction on k, this is 



9(a) 



«e^-i,z\r(m-l,n)\0 

But Pk-x,z\r{ m ~^-i n ) is isomorphic to i\z(w, ti) via the isomorphism a t— > aU{r}, 
and i a u{r] = i a + 1. Hence F k _ i z \ r (m — 1, n) equals (JTo|) . 

We can assume that the k — 1-assignment in Z is contained in the first k — 1 
rows and the first fc — 1 columns. Actually, we can assume that all elements of Z 
are contained in this fc — 1 by fc — 1 square, since if there is an element of Z in 
another column, say, then the row containing this element must be used in every 
fc — 1-covering of Z, and we can use induction on fc. 

We now subtract the minimum of the elements in columns fc, . . . , n, and sum 
over the m possibilities for the row r containing the new zero, using induction: 



F k ,z(m,n) = — + — V E 

m(n— k+1) m ^— ' 



3(a) 



i(n — k + 1) to ' (to — i a )(n — ja) 

r aGP fc ,z(m,n)\6 

Now we change the order of summation, and use the fact that the number of 
rows r for which r e a is i a . 

1 ^ i a g{a) 



F k , z {m,n) = — H V 



m(?i — fc + 1) 771 ^— ' (tti — i Q )(n — J a ) 

QGPfc,z(m,n)\0 

1 ^ 5(a) ^ g(a) 



E 3Vaj 



77i(?i-fc + l) ^ . (to — i a ){n - j a ) ^ .rrt(n-j a ), u 

aeP k: z(r,i,n)\0 aeP k , z (m,n)\0 (H) 

We have to show that the first term and the last sum cancel, that is, that 

3(a) = 1 

- (n — ja) 71 — fc + 1 
aePk,z(m.n)\0 

When a is the set of the first fc — 1 columns, we get a term which equals the 
right hand side. Hence we are going to show that everything else cancels. It clearly 
suffices to show that if J is a proper subset of the set of the first fc — 1 columns, 
then 

E 5(a) = 0. 

column— set (a) — J 

If J is nonempty, then the induced sub-poset of Pk t z(fn,n) consisting of those 
elements that contain J is isomorphic to Pk-\j\.z\,]( m , n— \J\) via the isomorphism 
a i— > a \ J. Hence the claim follows by induction on fc. 

If on the other hand J = 0, we have 

E 9(a) = E 9(a) - E 9(a). (12) 

aeP k ,z(m,n) a£P ktZ (m,n)\6 aeP k ,z (m,n)\6 

column- set(a)=0 column- set(a)^0 
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By what we already know, the second sum equals 1, since everything cancels 
except the term where a is the set of the first k — 1 columns. If we let cto be the 
intersection of all coverings of Z with k — 1 rows and columns, then 

9{a) = 9(a) = 1, 

a£P kz (m,n)\6 a<a 

by the definition of g. Hence 

J2 9(a) = 1-1 = 0, 

column — set(a)—0 

which proves the theorem. □ 



As a special case Theorem 3.7 enables us to compute the expected cost for 
completing a zero cost k— 1-assignment, thus answering a question posed in [ CS98| 



Corollary 3.8. Let Z k = {(1, 1), (2, 2), . . . , (k - 1, k - 1)}. Then 

(_l)fc-i-»-j 



F k ,z k (m, n)= Y 

i+j<k 



i,3,k- l-i-jj (m - i)(n - j) 



In particular, 



and 



k n — 7 



F k ,z k (k, k) = 1 (1 + 1/2 + • • • + 1/k) ~ 

Proof. The k— 1-coverings of Z k are all obtained by choosing, for each zero, whether 
to cover it with a row or with a column. It is easily seen that all subsets of k — 1- 
coverings (that is, all partial k — 1-coverings) can be obtained as intersections of 
fc— 1-coverings. The elements of the poset Pk,z k are therefore obtained by choosing, 
for every zero, whether to cover it with a row, with a column, or not to cover it. 
For any partial k — 1-coverings / a < /5, the interval [a, 0\ C Pk,z k is isomorphic 
to the Boolean lattice on |a| — \j3\ elements. We get 



/ i(6,i) = -^M«4)--E(- 1 ) H = 



Q>0 a>0 

fc-1 



, 1 )2' = -(1-2)^ = (-!)*. (13) 



1=0 



Any interval [0, a] C Pk,z k , a ^ is isomorphic to Pi,z n where / = k — \a\. 
Hence /i(6, a) — (— l) fc ~H , for all a € Pk,z k ■ The number of partial k — 1-coverings 
with i rows and j columns is given by the trinomial coefficient . fe _^^_,-). This 
proves the first statement. 



THE RANDOM ASSIGNMENT PROBLEM 



15 



For the second statement, observe that we can rewrite Fk i z h (m, n) by 

y*Y*( k - 1 \ {-if- 1 - 1 - 1 

j=0 i=Q 



i,j,k-l-i- jj (m - i)(n - j) 



Err hhn 



1 in — 7 ' — ' \ i I m — i 

j=o v J 7 J i=o v 7 

fe-i /fc-i\ 



[by Lemma O 1 = V / j /. • — -. (14) 

o (fc-i~-\) 



J=0 U--1-J 

In particular, when fc = m, we have 

1 fe_1 1 
i=o J 

□ 

4. Some consequences of the conjecture 
A consequence of the Main Conjecture is the following. 
Theorem 4.1. Let k <m,n be positive integers and Z the set of zero entries. If 



Conjecture l.i. is true, then there exist integers bij (independent of m and n) such 
that 

F k , z (m,n)= 7 ^7 v ( 15 ) 

im-i){n- 3 ) 

Moreover, if all the zeros in Z are in an ml x n' sub-array A! , then bi,j = 1 if 
i > ml or j > n'. 

Proof. Assume Conjecture [h^. Rewrite Fk,z{m,Ti) as 



E 



(m — x)(n — y) 

x+y<k V ;V U > 

where q x y is the probability that a random set with x rows and y columns (chosen 
uniformly) is a bad set. Now refine the counting of bad sets by letting s and t be 
the number of rows resp. columns that intersect A 1 — \l,m!\ x [l,n']. Then 

Ex T-^y (m — m'\(n—n'\ jk — l — y — x/ a/\ 
s=0 l^t=0 \ x-s I \ y-t ) a s,t \ A ) . 

" , (16) 



(m-x)(n-y) (™) (m - x) Q) (n - y) 

where d r st {A') is the number of sets of s rows and t columns that are not partial 
s + 1 + r- coverings. 

We use partial fractions with respect to m to get 

tm—m'\ mia{x,m' — 1} f-t \s—i (x\ ftn'— i+x— s— X\ 

\ 1 \ v V ) m f J (17) 



"x)( m - X ) fo 
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Inserting ( ft7| ) and the corresponding identity for A" f - into (16) we get that 



the contribution of 



(m— a;) (n — y) LlJ W K7 
a; 3/ 



to &ij is the integer 



(- i ) Hi (!)(l)EE(- 1 )^t'^:r 1 )te")<; 1 " ,, " I w. 

s=0 t=0 

if i < min{a;,m' — 1} and j < min{y, n/ — 1}. Otherwise the contribution is zero 
and the theorem follows. □ 



Theorem 4.1 enables us to write the value of Fk,z as a triangle of foil 's which 
has been very convenient and is the notation used in the Appendix. For example 



the resulting expected value in Example 3.3 can be written as 



-2 


1 


l| 


1 


1 




1 







The equation in the preceding proof can be used to give exact formulas for bij. 
Let r — x + y and sum over all x, y, then 

k — 1 m — k+r n — k+r / 

b id = i-{-i) i+ t J2 E E ((-i)' +t <7 1 - r (^ / )- 



r=i-\-j s—0 t—0 
r-max{j,t} 



r-max{i,s} 

jfe— 1— r / a t 



Here we have used the fact that d s t r {A') — Oii s + k— I— r> m' ovt + k— 1 — r > 



4.1. The probability that the smallest element in a row is used. In [092|, 
B. Olin conjectured, on the basis of experimental data, that in the case k = m = 
n, the probability that the smallest element in a row is used in the optimal n- 
assignment tends to 1/2, as n — > oo (Olin was considering uniform distribution in 
[0, 1] instead of exponential distribution, but this does not seem to matter in the 



limit). We now show that Conjecture 1.2 implies the conjecture of Olin. 



Note that the probability that the smallest element in a row is used is equal to 
the probability that the zero is used in a matrix with just one zero. This in turn 
is equal to the probability that the smallest element in a matrix without zeros is 



used. Assuming Conjecture 1.2, we can compute this probability. 



Theorem 4.2. Let A be an m x n-matrix with one zero. Let k be a positive integer. 



Assuming Conjecture l.i, the probability that the zero is used in the optimal k- 
assignment is 

1-®. 

mn 

Proof. We can assume that the zero is a^i, so Z = {(1, 1)}. Note that one element 
from the first row is always used in the optimal fc-assignment. We therefore compute 
the probability that another element in the first row, say ai^, is used. Let Z' = 
{(1, 1), (1, 2)}. Let pi t j and p- ^ be the corresponding probabilities for partial k— 1- 
coverings with i rows and j columns. We are interested in the difference pij —p'ij- 
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If i+j < k— 1, then pi : j = p' i • = 1. If i + j = k — 1, then Pij — p[ a is the probability 
that, if i rows and j columns are chosen randomly, position (1,1) is covered, but 
position (1, 2) is not. This is the case if and only if column 1 is used, but not row 
1 or column 2. The probability for this is 

to — i n — j j j(m — i)(n — j) 
to n n — 1 mn(n — 1) 

We have 



Pr[a h2 is used] = V - Pl ' 3 Pl ' 3 = V 

^ — ^ I m — lAim — i i — ' 



(m— i)(n— j) mn(n — 1) mn(n — 1) 

i+j<k y Jy J * i+j=k-l y ' y * 

Hence the probability that one of the elements ai^, ai,3, ■ ■ ■ , »i,n is used, is 

mn 

5. Some more evidence for the Main Conjecture 



□ 



If Conjectur e |l.2| is true then it should of course be consistent with the recursion 
from Theorem 2T in the case there is no doubly covered entry. In this section we 
will prove consistency of bij for the special case when i + j = k — 1. In this case 
the formula above simplifies to 

m'—l n' 

h^i-i-if- 1 y: £(- i r +t <*( A ') 

s=0 t=0 

One can of course do the same calculations starting with pij instead of qij and 
get that bij is 



m — s — 1\ i n — t — 1 



1 / \n' — j — 1, 
y instead of yiJ 



m'-ln'-l , , 1\ / ' + 1 

^■=(-i) fc - a e E(-w)fc;;: )(:: : , 

s =o t=o ^ m 1 L s \ n 3 L s (19) 

where g r s t (A') is the number of partial s + r + t-coverings of A' , with s rows and t 
columns. This is the formula we will use. 

We do the case when Z C [1,to'] x [1, n'] is optimally covered by n' columns. 
Let B x , x = l . . . ml be the matrix with zeros in ^ U (x, n' + l). Let also C be the 
matrix with zeros in Z U (to' + 1, n' + 1). Then the recursion step corresponding 
to the optimal covering of columns gives 

b .. {A) = E^liMgx) + (i-m')MC) = E^jAiCfM^Mgj | b (c) . 

' ' » ^(20) 

Note that we can include an extra row and column in A' if we just replace to' 
by to' + 1 and n' by n' + 1 in equation (|l^). To simplify the calculations we will 
do this and thereby assume that A' , B' x and C all are m' + l X ri' + 1 submatrices 
of A,B X and C respectively. 

We will now prove that equation ( |l9| ) remains valid during a step of the recursion. 
Partition g^X) = g°^(X) + g^(X) + g°J(X), by 

g®'t(X) = the number of partial s + i-coverings of X that use the last column. 
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g°' r (X) — ^ e num ^ )er °f P ar tial s + i-coverings of X that do not include the 
s ' last column but the last row. 

0,0/ y\ — the number of partial s + i-coverings of X that include neither the 
® s,t last row nor the last column. 

We can now establish a number of equalities from which the result will follow. 
First note g°J(C) = 0. Second g° s ^(A') + g%(A') = 9 %{C>) + g%{C>), this is 
used to match away the last term in (|2(]) and the "non-wanted" coverings of the 
enlarged A'. 

Third we have Y%Li 9°4{K) = m'g^iC')). 

We now put together 

T,7=i 9° s +i,t( B x) = Y^=i9 s°t( B 'x) = sg°s',t( A ')) and 
9s+i t t( C ') = 9 s ',t( A ')) to S et tne !ast identity needed 

Em (rn' — s\ 0,0/t->/\ (vn! — (s+l)\ / \-^m O.r / t-> \ / 0,r //~r\\ 

x=i\r n '-i)9..t\ B x)-{ m'-l ){2l x =i9s+i,A B x)-mg; +1:t {C)) 



0-0/ Al\ 

9 s ,t \ A )■ 



i 

' ml — 
m' — i / 

Plugging these identities and ( |l9| ) into recursion ( |20| ) everything will match 
nicely. 

6. Asymptotic consequences 
In this section, we show that certain limits can be computed from Conjecture 



1.2. We assume throughout this section that Conjecture 1.2 holds. 

Let if be a subset of the unit square [0, 1] x [0, 1]. Let n be a positive integer. 
Later, we will let n tend to infinity. 

We divide the unit square into squares Sij of side 1/n, such that 



1 i 



j ~ 1 j 



Let Z n be the set of all such that Sij intersects K. 

We shall show that under certain circumstances, we can compute the limit of 
F„ = F n> z n (n,n), asrn oo. 
For a fixed n we have 

Pij 



(n-i)(n-j) 

where pij is the probability that a set of i rows and j columns is a partial k — 1- 
covering. 

By the substitutions x — i/n, y — j/n, we can rewrite this as 

f 1 f 1 p n {x,y) , , 
Jo Jo (i-*)(i-y) dxd% 

where p n (x,y)=p y xn \ t y yn \ . 

If S is a subset of R 2 , we let M a ^{S) be the infimum of aX(Si)+bX(S2) (Lebesgue- 
measure), taken over all sets Si,^ of real numbers such that S C (Si xl)U 
(txS 2 ). 

Lemma 6.1. Suppose that K is compact. Let x and y be such that Mi— x ,i- y (K) < 
1 — x — y. Then p n (x 1 y) — > 1, as n — * oo. 
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Proof. We can find S 1: S 2 and <5 such that K C (Si x I) U (I x S 2 ), and (1 - 
x)X(Si) + (1 — y)X(S 2 ) < 1 — x — y — 8. By the compactness of K, we can assume 
that Si and S 2 are finite unions of intervals. 

Let n be a (large) integer. Then the number of intervals of the form [^-, ^] 
which intersect S± is at most n ■ A (Si) + constant, and similarly for the number of 
such intervals which intersect 5*2. Hence we can find a covering a of Z n using at 
most n ■ (A(S'i) + A(S , 2)) + constant rows and columns. 

Now choose a set a' consisting of \nx\ rows and \ny\ columns at random. The 
average number of rows and columns in a which are not in a' is at most 

n((l — x)X(Si) + (1 — y)X(S 2 )) + constant < n(l — x — y — 8) + constant. 

If n is large, then by the weak law of large numbers, the probability that \a — a'\ 
deviates from its mean value by more than n-S tends to zero. With high probability, 
|a(Ja/| < n. Hence with high probability, a 1 is a partial n — 1-covering. □ 

We say that a covering of Z n is minimal, if no proper subset is a covering of Z n . 

Lemma 6.2. Suppose that the number of minimal coverings of Z n grows slower 
than (1 + e) n for every e, as n — > oo. Let x and y be such that Mi_ Xt i_ y (K) > 
1 — x—y. Then p n (x, y) — > 0, as k — > oo. 

Proof. Let T(n) be the number of minimal coverings of Z n . Let a be a minimal 
covering. Since M\- Xt \- y {K) > 1 — x — y + 8, for a suitably chosen 6 > 0, we have 



for some Si slightly smaller than 8/2. 

Let a' be a random set constructed by letting each row belong to a' with prob- 
ability x, and each column with probability y. Since with probability at least 1/4, 
a' will contain at least i rows and j columns, in order to prove that p n (x,y) tends 
to 0, it will suffice to show that the probability that one can cover Z n by adding 
n — i — j — 1 rows and columns to a' tends to 0. 

Let B(p) denote the random variable that is 1 with probability p and with 
probability 1 —p. The number of rows in a\ a' is a sum of i a independent B(l — re- 
distributed random variables, and similarly, the number of columns in a \ a' is a 
sum of j a independent B(l — j/)-variables. In order that the size of a \ a' be at 
most n — i — j — 1, one of these sums must deviate by at least n ■ Si from its mean 
value. By standard estimates in probability theory, this is at most e~ c ' n , for some 
constant c depending only on 8±. 

Hence the probability that one can obtain a covering of Z n by adding at most 
n — i—j — 1 rows and columns to a' is at most T(n) ■ e~ c ' n , which, by assumption, 
tends to as n — > oo. □ 

Corollary 6.3. If K is compact, and the number of minimal coverings of Z n grows 
slower than (1 + e) n for every e, as n — > oo, then 

f dx dy 



(1 - x)i a + (1 - y)j a > n ■ (1 - x - y + 8). 
If we let i = [n ■ x\ and j = [n ■ y\ , we certainly have 

(1 - x) ■ i a + (1 - y) ■ j a > n - i - j + 2Si ■ n, 




where D is the region in which Mi_ x ^_ y (K) < 1 — x — y. 



20 



SVANTE LINUSSON AND JOHAN WASTLUND 



6.1. An example: Zeros outside a quarter of a circle. Let K be the region 
x 2 + y 2 > 1, < x, y < 1. To find M a ^(K), for given a and b, we have to find the 
point where the slope of the curve x 2 + y 2 = 1 is —a/b. The part of the curve which 
is above this point should be covered with a horizontal strip, and the part which is 
to the right of this point should be covered with a vertical strip. This point is 

a b 



yja 2 + b 2 y/a 2 + b 2 



Hence 



M a , b (K) = a 1 - 



Va 2 + b 2 



6 1- 



Va 2 + b 2 



which simplifies to 



a + b - \/a 2 + b 2 . 

The number T{n) of minimal coverings of Z n is at most n+ 1, since any minimal 
covering consists of, for some i, < i < n, columns i,...,n, and all rows that 
contain some element of Z n not covered by the columns. 

The region D is given by the inequality M{K\- x ^i- y ) < 1 — x — y, which becomes 

(1 - xf + (1 - y) 2 > 1. 

In other words, D is the part of the unit square which is outside the circle of radius 
1 centered in the point (1,1). By Corollary 6.3, we have 

dxdy 



lim F n = 



D 



(l-x)(l-y) 



71" 

24' 



6.2. Generalizing the example. Fix p > 1. Let K be the region given by 
x p + y p > 1, < x, y < 1. We wish to compute M a ^(K). To find the point on the 
curve x p + y p — 1 where the slope is —a/6, we compute: 

^ = -a*- 1 ^ - x^ 1 ^- 1 = ~{x/y) p -\ 
ax 



It is not difficult to see that the point is 
/ 



Of" 1 



bp- 1 



( _E_ ,_E_\ 1 / p ' / _E_ P \ 

Map- 1 + op- 1 J la?- 1 +0P- 1 ) 



i/p 



We have 

M a , b (if) = a 



1 - 



(Xjj-i 4- 6 p — 1 
p 

or- 1 



fop- 1 



i/p 



b p— 1 



/ — E 
I Of" 



1 + bp- 1 



i/ P 



clp- 1 + bp- 1 



i/p 



a p-i -)_ p-i 



(21) 



Hence the inequality Mi- x _i^ y (K) < 1 — x — y becomes 



((1 - x)p- x + (1 - y)p- 1 ^j 



> 1. 
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Again the number of minimal coverings of Z n is at most n + 1. To find the limit of 



F n , we now by Corollary 6.3 wish to compute the integral 

' 1 dx dy 



c =o Jj/=(i-a;") 1 /" xy 
where u — p/(jj — 1). We eliminate the inner integral, and obtain 
i-log((l-^A^ = _l ^log(l-^) 



/ *^ ^ J ^ 

We make the substitution t — x u , which gives dx — — ^=r. We get 



1 /Mo g (l-t) dt= _l [ l (-t-*/2-*/B-. = 



t u 2 J \ t 



1/11 \ 1 7T 2 , , 

1 + 7 + 7T + --- =7^-^T- (22) 



4 9/ u 2 6 



Substituting back pj (p — 1) for u, 

lim F„ = I 1 - 



n — >cx: 



1 \ 7T 



p J 6 



7. Remarks 

Looking at the triangles of frjj's one discover certain patterns. Define the prop- 
erty acyclic recursively by first letting the empty set be acyclic. Second, let Z be 
an acyclic set that has an optimal covering a with only rows or only columns. Let 
be any position not covered by a. Then we define ZU to also be acyclic. 



Corollary 7.1. Assume Conjecture 1.2. If Z is acyclic, then 6 ,o = ( — 1)'^' ( \z\)- 

□ 

This corollary can be proved by first specializing (18) to &o,o an d then study 
the set of bad sets and find the proper bijections which proves it by induction on 
\Z\. The proof we have found is much in the spirit of the proof in Section |^ of 
consistency for b^j, i + j = k — 1. There are some highly non-trivial details and we 
omit the proof. 

Given a set of zeros Z, let Z row (i), be the number of zeros in row i and Z co i(j) 
the number of zeros in column j. Ordering the nonzero Z row (i), 1 < i < m by size, 
we get an integer partition of \Z\, call it X r0 w(Z). Similarly we define X co i{Z). 



Conjecture 7.2. Given k and Z, write ffc,z(w, n) = X^+j<fc (ro-i)'(rt-j) f or some 
(uniquely defined) integers &$j as in Theorem ^. \ . Then the value o/&o,0; &i.o> ■ 
only depend on \ row (Z) and bo,o, &o,ii • • • j &o,fe-l only depend on X co i(Z). 



'fe-1,0 



It is not difficult to see that this conjecture implies Corollary 7.1 



We mention that Theorem 2.10 can be generalized in the following way: 



Theorem 7.3. Let at and bi, 1 < i < N, be real numbers. Let X be an ex- 
ponentially distributed random variable with mean 1. Define the random variable 
L G {1, . . . , N} by ai + b[X = mhij ai + biX . Then E(bi) = E(aj + b[X) — min(ai). 
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Proof. Let f{x) = mini(ai + bix). Then 

f°° d 

- mm( ai ) = -/(0) = jf — (e- x f(x)) dx = 



e' x f{x)dx- I e~ x f(x)dx = 
Jo 

= E(/'pO) - E(/(X)) = E(6j) - E(o/ + bjX). (23) 

□ 



As a last remark we want to point at the recently published preprint [BCR00|, 
where they generalize the random assignment problem in a different direction. They 
assume that the matrix has independent exponential random variables with differ- 
ent intensities. In the case when the intensities form a rank 1 matrix they give 
a conjecture of the expected value of the optimal ^-assignment. Observing that 
the binomial coefficient in their conjecture is the Mobius function of a truncated 
Boolean algebra, it is not difficult to join our conjecture with theirs to the case 
with a rank 1 matrix of intensities, but in some entries the random variables are 
replaced by zeros. We have done no computations to see if this join of conjectures 
is true. 
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Appendix: Computer generated verification of Conjecture 1.1 

k < 5 



FOR 



Our algorithm for calculating Fk y z(m,n) is essentially a formalization of the 
method outlined in Section |[ It was implemented in Maple to take advantage of 
Maple's built in functions for simplifying rational expressions. 

Let us immediately state that we do not have a proof that this algorithm always 
returns a result. However, if a result is returned, we know that this result is correct. 
The difficulty lies in that this algorithm does not use fixed values for m and n, but 
tries to compute Fk(m,n) as a rational function of m and n. For fixed numerical 
values of m, n and k, it is in principle always possible to calculate F^(m, n). 

The algorithm uses a data-structure which we will here call an array. An array 
is a matrix whose elements are finite sums of the form 

i 

where X\, X%, X3, . . . are independent exp(l)-distributed random variables, and 
is a rational function in the variables m and n. An array M is thought of as 
representing a window of an m by n matrix A, all of whose entries outside the region 
visible in M are exp(l)-variables, independent among themselves, and independent 
of all entries in M. In this way we can for example represent a pattern Z of zeros, 
independently of m and n, and compute -Ffc,z(m, n) as a rational function in m and 
n. 



Below are the cases computed to calculate F^(m,n) using Theorem 2.1. The 
same computer program has been used to verify Conjecture [0] also for k = 6, but 
this computation needed about 10,000 cases, and is not written out here. In each 
case, M is an array, and an optimal covering of the zeros of M is computed. The 
other cases that are used for a particular case are listed. The expected value of 
the optimal ^-assignment is denoted F. The value of F is given as a matrix. This 
matrix represents the coefficients bij. Hence the entry in position (i + X,j + 1) of 
the matrix is the coefficient of 1/ ((m — i)(n — j)) in the expression for F. Sometimes 
the value of F cannot be written as a linear combination of such terms. Then the 
extra terms are written out explicitly. 

8.1. The case k = m = n = 7. In case k = rn = n, we can take advantage of the 
fact that every fc-assignment contains one element from every row and every column. 
This means that we can start by subtracting the minimal element from every row, 
and then continue by subtracting the minimal element from every column that does 
not already contain a zero. Hence we need only take into account the cases where 
Z contains at least one element in every row and every column. It turns out that 
in the calculation of Fr(7, 7), all the cases that occur can easily be reduced to cases 
that have already occurred in the calculation of Fk(m, n) for fc < 5. In this way we 
have verified that F 7 (7, 7) = 1 + 1/4+ 1/9 H h 1/49, as conjectured. 
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Case 1 

k = 0, M = [], F = 



k = 1, M 



using case 1 



covering row 1 
using cases 1 2 



k = 2, M = [ ] , F = [ J I 



k = 2, M 



using case 3 



[ X! ] 



1 




covering rows 1 2 
using cases 12 3 



111 

-2 10 
10 



covering column 
using cases 3 5 



covering rows 1 2 
using cases 1 3 



3, M = ^ j , F = 



using case 8 



covering rows 12 3 
using cases 12 3 5 









x 2 
x 4 



> F -- 



covering row 1 
covering column 1 
using cases 3 5 10 








*1 







k = 4, M = 








^2 









X 3 


X 4 _ 





covering rows 12 3 
using cases 13 7 



Case 13 
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covering columns 1 2 
using cases 7 11 12 



covering rows 12 3 
using cases 12 5 



covering columns 1 2 
using cases 3 10 14 




*4 



[S 







covering rows 1 2 
using cases 13 15 16 



: 4, M : 





x 2 



covering columns 1 2 
using cases 7 11 16 



fc = 4, M 



-[S] 



covering column 1 
using cases 17 18 



covering rows 12 3 
using cases 1 7 



x 2 







x 2 



covering rows 1 2 
using cases 18 20 



covering row 1 
using cases 19 21 



Case 23 
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(m — 2) n 



X 2 
X E 





X-L 
x 3 

x R 



(TO - 2) (2 m n - 5 — 1) 



covering rows 12 3 
using cases 12 3 5 



■ 2)(n - 1) 
























— 1 








Xl 










-1 


1 





*3 




, F = 




-1 


— 1 
1 


1 








*5 


*6 


*7 - 




1 












2 (m 2 - 4m + 4)(2mn-2m- 



. +4) 



2 (m 2 — A m -\- 4) n 



covering rows 12 3 4 
using cases 1 2 3 10 24 
























Xl 








x 2 


x 3 





x 4 


x 5 


x 6 . 



covering rows 12 3 4 
using cases 1 2 3 5 10 



: 5, M : 



covering rows 12 3 4 
using cases 1 2 3 5 24 



Xj 

2) (n - 
2 m - 




X 2 



Xl 

x 3 



x 5 x q 



2 (m 2 -4m + 4)(2mn-2m- 



1 +4) 



2 (m 2 — 4 ?7i + 4) n 



3, JVf : 



1) n 




*1 



1 




(to — 1) (2 m n - 



covering rows 1 2 
using cases 12 3 




X! 



1) („ 





*2 



2 („ 



2 m + 1) (2 m 7r - 



. + 2) 2 (m 2 - 2m + 1) 1 



covering rows 12 3 
using cases 28 1 2 5 



X 2 



1) (n - 1) 




Xl 

x 3 



2 („ 



> + 1) . 



2 (m 2 - 2 to + 1) (2 r 



• + 2) 
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covering rows 12 3 
using cases 28 1 2 3 



: 5, M : 






A'3 



(m - 



2) („ 










o - 










1 


-2 


1 







x 2 







1 


1 


1 








x 4 


*5 


, F = 


-1 
-1 


-2 
1 


1 











x 6 


x 7 


x 8 . 




1 















> + 4)(m-l)(2r 



i + 4) (m 2 



covering row 1 

covering columns 1 2 

using cases 29 30 3 5 10 25 26 27 

Case 32 














o - 













— 1 


1 












x 2 







1 


— 2 


1 







-^3 





X 4 


*5 


, F = 





-2 


2 










(777 - 2)(„ - 1) 








-1 


1 
















*6 


x 7 


*8 - 




1 














2 777 - 


- 5 






1 


1 










( m 2 _ 4 771 + 4) (2 m 


77 - 2 7T 


- 5 r 


+ 4) 


2 (7772 


- 4 777 


+ 4)77 









covering rows 12 3 4 
using cases 1 2 5 14 24 









x 2 




x 1 o - 

x 3 
x 4 x 5 _ 










(777 - 


- 1) (« 




- 1) 


, F = 










1 

-2 


— 2 

2 


1 




1 




1 


1 


— 1 
1 


1 











+ 2 (ml - 2 


777+l)(2 77777-27T 


- 377 + 2) 2 ( 


777 2 - 2 771 4- 1) 77 



covering rows 12 3 
using cases 28 1 3 7 

Case 34 



: 5, M : 























Xl 
















- 1) 





x 2 


x 3 




(' 


,, - 2) (71 - 












x 4 







*5 


*6 











x 7 















Xio 


Xn 


x 12 - 







1 


-2 


1 - 









1 1 


1 














-2 1 








+ 2 







1 










171 — 2771 — 571 + 4) (777 



















covering row 1 
covering columns 1 2 
using cases 28 3 5 10 24 

Case 35 





















Xl 


x 2 





*3 


x 4 


X;, 





*6 


x 7 


*8 



covering rows 1 2 
covering column 1 
using cases 31 32 33 34 7 14 

Case 36 



: 5, M : 









Xl 














x 2 








x 3 


x 4 





*5 


x 6 


x 7 



covering rows 12 3 4 
using cases 1 3 7 12 



Case 37 
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*1 


a 











x 2 







x 4 


x 5 





x 6 


x 7 





covering rows 1 2 
covering column 1 
using cases 36 3 5 10 12 26 



covering rows 12 3 4 
using cases 1 3 5 7 16 

























x 3 





x 4 





covering columns 12 3 
using cases 35 37 38 16 



11 




















*1 














*5 


*6 


x 7 _ 



covering rows 12 3 4 
using cases 1 2 3 10 























*2 


x 3 





x 4 


*5 





*6 


*7 


*8 



covering row 1 
covering columns 1 2 
using cases 40 3 5 10 14 26 



1) (n + 1) 



■ 1) (2 m 7 



-2 + m) 



1) n 



X 4 



> + 1) (2 r 



covering rows 12 3 
using cases 42 1 3 7 





x 2 



■ 1) (n - 


x 7 







. + 1) (n 



- 2 m + 1) (2 » 



- 3 n + 4) 



covering rows 12 3 4 
using cases 42 43 1 3 12 



Case 45 
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o 

x 4 

■ 1) (n - 1) 

2 m - 3 



(m 2 - 2 m + 1) (2 r 

covering rows 12 3 4 
using cases 42 43 1 3 7 



*1 
*2 



■ + 4) 



> + !)(« 



(771 - 1) (77 - 1) 





x s 



x 2 
x 4 

Xfi 



covering columns 12 3 
using cases 36 41 44 45 7 



- 2 m + 1) (2 m t 



in + 4) 



2 m + 1) (n - 1) 















11 





^1 


x 2 


*3 


x 4 





*5 


11 




x 7 


*8 



-3 
2 




covering rows 12 3 4 
using cases 1 2 5 10 14 



covering rows 12 3 
using cases 42 1 2 5 

























*3 




x 4 


(rr 


-1)17 










-1 


1 









-1 


1 









-1 


1 










+ 


1 
















2m + 1) (27 



- + 1) 



(m2 



. + 1)7 



5, M 













x 4 

*6 





*1 



1) („ 


x 7 




*5 



covering rows 12 3 4 
using cases 42 48 1 2 10 



- 2 m + 1) (2 to. 7 



- 3 n + 4) 



+ 1) (n 


















*1 




*2 







*3 



( m 2 



■■+!)» 



2m + 1) (2 7 



■377 + 1) 



covering rows 12 3 
using cases 42 1 2 3 



X 2 
1) („ 







- 2 777 + 1) (2 777 7 



3 777 - 3 77 + 4) 



: + 1) (77 
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covering rows 12 3 4 
using cases 42 50 1 2 3 



covering columns 12 3 
using cases 47 49 51 3 26 



1) (n - 1) 



X 4 



*2 




I) 2 (r. 



>(m - l) 2 (2, 



■ + 4) 



*5 

x 7 



J rows 12 3 4 
using cases 1 3 7 10 16 



k = 5, M : 



(m - 1) (n - 1) 




x 2 
X3 X4 



( m 2 - 2m + l)(« 



- 1) 
(2 m- 3) 2 



(m 2 - 2 m + 1) m (3 m - 4) (2 m n - 



- J» + 4) 



covering rows 12 3 
using cases 39 46 52 53 



k = 5, M = 




*4 

2) („ 




covering rows 12 3 4 
using cases 1 3 7 16 24 



2 (m 2 - 4m+4)(2mn-2m- 



1 +4) 



2 (m 2 — 4 ?7i + 4) n 






*1 














*2 


*3 







*5 


x 6 





*7 


*8 


*9 



iws 1 2 
tlumn 1 

1 31 33 34 55 7 16 



5, M 



1) (n 


*5 



x 10 



> + l)(n 



. + 1) (2, 



: + 4) 



covering rows 12 3 4 
using cases 42 43 1 5 16 
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i) („ 
o 





x 5 



A'3 
*6 



(2 m- 3)^ 



X 4 

x 7 



(3 to - 4) (to - I) 2 (2 to n - 3 m - 3 n + 4) 



(3 to - 4) , 



vs 1 2 3 

41 56 57 16 



: 4, JVf : 



*5 



X 6 X 7 



(6 - 5 to + i 



2)(2, 



covering rows 12 3 
using cases 1 2 3 5 24 



covering rows 12 3 4 
using cases 59 1 2 3 10 



X 2 + ■ 



+ - ■ 

2 

3 



■ 2) (n - 1) 
) 

(2 m - 5)" 









o - 




" 








— 1 


1 








X-L 










-1 


1 








x 4 


*5 


, F = 






-1 
1 


1 












x 6 


x 7 


*8 - 




1 















12) (2, 



2 („ 



+ 4) 



• + 4), 



5, M -- 



X 4 + ■ 



(to - 2) (n 



X 2 
*6 




*1 

x 3 
x 7 



2 („ 



2 („ 



• + 4), 



2 (m 2 



. + 4) (r, 



■3) (2m. 



■ + 4) 



covering rows 12 3 4 
using cases 59 1 2 3 5 



; 3, M ; 



X 2 + ■ 



*3 




Xl 



(2 TO T 



!)(„ 



1) (m - 2) 



covering rows 1 2 
using cases 28 1 2 3 



: 4, M : 



*1 + 



*2 



(to - 1) („ 



X 4 





x s J 



, + 1) (r, 



2) (2, 



• + 2) 









-1 1 





-1 


1 





1 





1 








1 




- 2 


2 ( 


™2 _ 2 


m + l)n 



2 („ 



covering rows 12 3 
using cases 62 1 2 5 
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: 4, M - 



covering rows 12 3 
using cases 62 1 2 3 



2 („ 



X;, 



- 1) (" 
1 

+ ■ 





x, 



x 4 



» + 1) n 2 (m - 2) (n 



- 2 m + 1) (2 7 



■2 m - 3 ti + 2) 



2 („ 



X 3 + • 



x 4 



(m - 2) (n - 1) 










o - 




a 





1 


-2 


1 - 







Xl 


X 2 







i 


1 


1 





3 





X 5 


x e 


, F = 


-l 


-2 


1 















1 











2 


x 7 


X 8 


x 9 - 




l 

















2 (m-l)n (m^-4m + 4)n 
{1m. ~ 5) 2 



(m 2 - 4 m + 4) (m - 1) (m - 3) (2 m n - 2 m - 5 n + 4) 



covering row 1 

covering columns 1 2 

using cases 60 61 63 64 3 5 10 26 



X 3 + 



X 4 



(m - 2) (n 










o - 













— 1 


1 


Xl 





X 2 







1 


-2 


1 








X 5 


X 6 


, F = 






-2 
1 


2 











x 7 


x 8 


X 9 . 




1 















2 (16 ? 



■ 12) (2, 



. + 4) 2 (r. 



• + 4), 



2 (m — 3) n 



covering rows 12 3 4 
using cases 59 1 2 5 14 



X 2 + • 



X 3 



(m - 1) ( 




Xl 



1) 

(2 m- 3) 2 





x 4 

x 6 



, + 1) (2, 






1 


-2 1 





— 2 


2 





1 





1 








1 




n - 2 


2 ( 


Tl 2 - 


2 m + 1) n 



covering rows 12 3 
using cases 62 1 3 7 



: 5, M : 



^1 + 



covering row 1 
covering columns 1 2 
using cases 59 62 3 5 10 





Case 68 
















X 2 











x 3 


X 4 


(m _ 2 ) („ - 1) 








X 5 





x 6 


x 7 





X 8 


x 9 


x 10 





Xll 


x 12 


x 13 



■ 2) (n 



1) (>7i - 3) (2 r 



. + 4) 












a 


Xl 





x 2 


x 3 





x 4 


x 5 


x 6 





x 7 


Xs 


x 9 _ 



— 2 
1 






covering rows 1 2 
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covering column 1 

using cases 65 66 67 68 7 14 














X 3 


, F 


X A 


X 5 




x e 


X 7 - 





■ 2 m + 1) (n 



covering columns 12 3 
using cases 37 57 69 16 



i + l) (2. 



*4 + ■ 



2) („ 



. + 4) 



Cast 


71 







Xl 








*2 




X 6 





x 7 


*8 


*9 


x 10 



2 (771 - 3) 1 



2 („ 



1 + 4). 



2 (77 



>+4)(t> 



■ 3) (2 , 



. + 4) 



covering rows 12 3 4 
using cases 59 1 3 7 16 






Xl 








*2 





X 3 


X 4 





X 5 


X 6 


x 7 





^8 


X 9 


x 10 



iws 1 2 
>lumn 1 

; 65 67 68 71 7 16 






a 







3 


1 


1 


1 










X 2 


, F = 




1 


1 











x 3 


X 4 _ 




1 


1 



















1 















covering row 1 
covering column 1 
using cases 54 58 70 72 






*1 


X 2 








x 3 





x 4 










x 6 





x 7 


x 8 


x 9 . 



covering rows 12 3 4 
using cases 1 7 20 





X 2 



covering columns 12 3 
using cases 56 74 20 



covering rows 12 3 4 
using cases 1 3 5 12 16 
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x 2 





X- 


x 4 





x 5 


x[ 





x 7 


x 8 


x < 







covering columns 12 3 
using cases 36 38 41 76 7 
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Xi 





x 2 


o x 3 













x e 


X 7 Xg _ 



covering rows 12 3 4 
using cases 1 3 5 7 16 



5, M . 




X 3 




Xl 


x., 



covering columns 12 3 
using cases 36 41 76 78 7 



: 5, M ■■ 






Xl 


X 2 













x 4 


x 5 





x e 


x 7 





X 8 


Xg 


x 10 



ws 1 2 3 4 
1 7 12 20 



X 4 



covering rows 12 3 
using cases 75 77 79 80 



: 5, M 






*1 





x 2 





*3 





x 4 


*5 







X 7 _ 



• lumns 12 3 
56 72 80 20 





x, 



covering columns 1 2 
using cases 73 81 82 



covering columns 12 3 
using cases 37 69 76 16 



k = 5, M = 















Xl 


x 2 





X' 


x 4 





xa 


x\ 
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